Test FAMS on 5D data
Contents
randomize data
vect1 = rand(1000,5)*100;
sigma2=[1,3,1,3,2]*5;
mu2=[70,60,70,55,89];
vect2 = randn(120,5);
vect2 = vect2.*repmat(sigma2,size(vect2,1),1)+repmat(mu2,size(vect2,1),1);
sigma3=[1,3,1,3,2]*5;
mu3=[20,10,15,55,9];
vect3 = randn(100,5);
vect3 = vect3.*repmat(sigma3,size(vect3,1),1)+repmat(mu3,size(vect3,1),1);
data=[vect1; vect2; vect3];
run FAMS
[modes, p_modes, w_modes]=fams(data,20,15,200,'5Ddata','res/','-h',20);
Load data points from matlab ...done
RunFAMS with res/pilot_200_5Ddata.txt ...
Running FAMS with K=20 L=15
Run pilot fixed bandwith...done.
Start MS iterations..........done.
Join Modes with adaptive h/2, min pt=40, jump=1
pass 1.done
pass 2..........nrel 6
0 0
1 1
2 2
3 3
4 4
5 5
done
Save convergence points ...done
Save joined convergence points ...done
Save indicies of modes ...done
FAMS done.
view results
format short g;
p_modes
p_modes =
446 70.369 57.104 69.913 55.933 88.223
315 21.429 8.8381 14.123 56.03 8.5525
98 76.775 26.421 71.663 52.89 48.564
58 84.778 69.555 45.305 78.922 29.83
54 44.767 34.338 36.552 46.509 30.843
46 27.184 33.584 54.598 27.461 33.91